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ABSTRACT 

A solution to a given equation is structurally stable if it suffers only an infinitesi- 
mal change when the equation (not the solution) is perturbed infinitesimally. We 
have found that structural stability can be used as a velocity selection principle 
for propagating fronts. We give examples, using numerical and renormalization 
group methods. 
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The steady state equation for a travelling wave propagating into an unstable 
state does not always uniquely determine the wave speed. Instead there may be 
multiple stable steady state travelling wave solutions, even though the physical 
system described by the equation exhibits reproducibly observable behavior cor- 
responding to only one of these solutions. 1 6 In such a situation, it is desirable to 
formulate a so-called selection principle, which would allow one a priori to distin- 
guish observable from unobservable steady state front solutions without having 
to solve directly the equation of motion starting from the initial conditions. 

For a certain class of equations, rigorous analysis shows how a wide range 
of physically realizable initial conditions evolve into the selected front, which 
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turns out to be the slowest stable solution allowed by the steady state equation. 
A physical, heuristic interpretation of this result, known as the linear marginal 
stability hypothesis, has been proposed and is believed to be applicable in the 
so-called pulled case, for which the selected speed may be determined by the 
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linear order terms alone. ' However, it is well-known that there is another case, 
the so-called pushed case, where analysis of the linear order terms alone is not 
sufficient to determine the speed, and the linear marginal stability hypothesis 
fails. 8 ' 9 

The purpose of this Letter is three-fold. First, we recall the notion of struc- 
tural stability — the stability of a front with respect to a perturbation of the 
governing equation — and argue that only structurally stable fronts are observ- 
able. We next show that for structurally stable fronts, a renormalization group 
(RG) method can be used to compute the change in the front speed when the 
governing equation is perturbed by a marginal operator. Finally, by combining 
the structural stability principle with RG, we are able to predict the selected 
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front itself. Our results apply to both the pulled and pushed cases. Roughly 
speaking, structural stability is an insensitivity to model modifications, whereas 
the RG may be interpreted as a method to extract the structurally stable be- 
havior of a model. 10 ' 11 Structural modifications of travelling wave equations have 

12 

been studied previously (e.g., ZePdovich's work on flame propagation ) but to 
our knowledge, structural stability has not previously been proposed as a selec- 
tion mechanism. RG methods have previously been used to study the asymp- 
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totics of partial differential equations (PDEs) ' and propagating fronts in the 
Ginzburg-Landau equation. 14 

A good model of reproducibly observable physical phenomena must give 
structurally stable predictions. That is, the observable predictions provided by 
the model must be stable against "physically small" modifications of the system 
being modelled. We will quantify below the meaning of the term "physically 
small" for a certain class of reaction-diffusion systems. The idea of structural 
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stability used here is close to that proposed by Andronov and Pontrjagin for 
dynamical systems. In the modeling of natural phenomena, we need not require, 
as did Andronov and Pontrjagin, the structural stability of the entire model, but 
need only to require it of the solutions corresponding to reproducibly observable 
phenomena. We call these structurally stable solutions. Our structural stability 
hypothesis states that only structurally stable solutions of a model represent re- 
producibly observable phenomena of the system being modeled. This hypothesis 
is implicit in most mathematical modeling, and indeed often redundant, yet we 
will demonstrate that for reaction-diffusion equations, this hypothesis correctly 
singles out observable propagating fronts. The basic reason for its efficacy in 
the situations studied here is that the formulation of reaction-diffusion models 
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sometimes inadvertently includes an unphysical feature, although the model is in 
some sense close to a class of physically correct models. 

Consider Fisher's equation 1 on the interval — oo < x < oo: 

dip d 2 ip 

where F is a continuous function with F(0) = F(l) = 0. We will usually be 
interested in boundary conditions where ip is zero at one boundary and unity 
at the other. If F satisfies the condition: F{ip) > for all ip G (0,1), then 
there exists a stable travelling wave solution interpolating between ip = 1 and 
ip = with propagation speed c for each value of c greater than or equal to some 
minimum value c*. The positivity condition on F stated above together with 
differentiability of F at the origin will henceforth be called the AW-condition; 
when it is satisfied, c* > c = 2y / F'(0) . Aronson and Weinberger 7 proved that 
for (1) with the AW-condition satisfied, the selected solution is that with speed 
c*. In most systems studied by physicists, the minimum wave speed satisfies 
c* = c, which corresponds to the pulled case. Often, the initial conditions decay 
sufficiently fast (faster than some exponential function) to ip = that the selected 
wave speed is in fact c*. The pushed case is equivalent to the statement c* > c. 
In this paper, we are concerned not only with Fisher's equation subject to the 
AW condition, but with other equations or systems of equations not satisfying 
the conditions required for Aronson and Weinberger's rigorous proof, but which 
still exhibit the selection problem. 

It is straightforward to show that all propagating solutions of (1) are struc- 
turally stable against C 1 -small perturbations 5F of F. Unfortunately, reaction- 
diffusion equations are not in general structurally stable with respect to C°-small 



perturbations. Consider (1) as describing the propagation of fire along a fuse. F 
represents the net rate of heat production as a function of temperature ip. The 
value ip = corresponds to the flash point, and ip = 1 corresponds to the steady 
burning temperature. It is reasonable that the observable properties of such a 
front would be insensitive to most small changes to F. However, by altering F 
very near ip = with a C°-small perturbation, dF/dip in the neighborhood of 
tfj = can be made arbitrarily large. That is, the rate at which heat production 
increases as a function of temperature at or near the flash point can be made 
very large, and this explosive low temperature behavior will travel very rapidly 
along the fuse. 

It is clear then that certain C°-small perturbations are not physically 
small. This is the case, however, only for perturbations which increase 
su -PtpG(o,ri}(F(ip) /ip) appreciably for some r\ > 0. We will call a C°-small perturba- 
tion for which sup ^ >0 (8F (if;) /ip) is less than some small positive number (which 
goes to zero continuously as the C°-norm of SF vanishes) a p-small perturbation. 16 
The precise form of our structural stability hypothesis is: physically realizable 
solutions of (1) are those which are stable with respect to p-small structural 
perturbations. 

The ordinary differential equation (ODE) governing the travelling wave front 
shape = ip(x, t) can be transformed into the equation 

dU 

v =- cv -— (2) 

with the identifications £ = x — ct — > t, ip — > q, difj/dt; — > q = p, and F = dU/dq. 
This ODE describes the position q of a unit mass particle subject to a potential 
U (q) and friction. The coefficient of friction is c, the speed of the travelling wave. 
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Traveling- wave solutions of (1) interpolating between the fixed points ip = 
and ip = 1 correspond in this particle analogy to trajectories which begin at the 
maximum of U located at q = 1, with zero kinetic energy, and which terminate at 
the origin. Those trajectories which correspond to stable solutions of the original 
PDE are those for which q never changes sign. 

If the origin is not an isolated local minimum of U, there is only one value of 
the coefficient of friction which allows the particle to stop here without overshoot- 
ing. If the origin is an isolated local minimum of the potential, there is a critical 
value c* of the frictional coefficient. For all smaller values, the particle overshoots 
the origin at least once, while for a continuous set of larger values, it converges 
to the origin as t — > oo without overshooting. For both types of systems, we 
define c* to be the smallest value (unique value in the former case) for which the 
particle approaches the origin in the t — > oo limit without overshooting. 

We have proven that c* is continuously dependent on the continuous func- 
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tion F. Thus because we can make the origin a local maximum of U with 
arbitrarily small C°-perturbations (in fact, p-small perturbations), and because 
c* is the only value of c for which the particle stops at the origin (without over- 
shooting) in this case, the continuity of c* implies that only this critical value is 
structurally stable. Therefore, our structural stability hypothesis asserts that c* 
is the unique observable front propagation speed of the original PDE. 

For AW-type equations, it has been proven that c* is indeed the selected 
speed. Our structural stability hypothesis is thus correct in this case. For 
non AW-type equations, however, there is no proof that this minimum value 
is selected. We have therefore performed numerical studies to test our conclu- 
sions for non-AW equations as well as for systems of coupled reaction-diffusion 
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equations. 17 ' 18 An example is the equation studied by van Saarloos: 20 



where 7 < 1/12. For < b < 1/2, the front is pushed, whereas for 1/2 < b < 1 the 
front is pulled. We replaced the potential term by 6(ip—A)(ip—A)(l—ip)(ip+b)/b, 
where 9 is the step function, and let A — > + . For both pulled and pushed cases, 
the unique selected velocity converged to that dynamically selected. 

We now proceed to show how RG can be used to compute the change in the 
front speed when an equation, whose structurally-stable exact solution is known, 
is perturbed by ap-small operator. Introducing new variables X = e x and T = e*, 
the propagating front solution reads ip(x — ct) = $(XT~ C ). Thus the front speed 
is interpreted as an anomalous dimension. The renormalization group theory 
applied to PDEs 11,13,1 should therefore be applicable here, too. In terms of RG, 
what we have found above is: p-small C°-perturbations are at worst marginal 
perturbations, but generally, C°-perturbations can be relevant (actually, in some 
sense much worse, since the changes they produce are in certain cases indefinitely 
large) . 

Let ipo(x — cot + xq) be a stable travelling front solution of (1) with speed 

Co and constant of integration xq. Let us add a p-small structural perturbation 
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5F to (1), where its sup-norm \\SF\\ is of order e, a small positive number, 
and assume that in response the front solution is modified to ifjo + Sip. Defining 
£0 = x — cot + xo and linearizing (1) with respect to e in the moving frame with 

velocity cq, we formally obtain the following naive perturbation result: 

t +00 

<ty(£o,f) = e-^y dt> J di'G^^ty^'lHF^'))- (4) 

to —00 



Here to is a certain time before 6F(ipo(£o)) becomes nonzero, and G is the Green's 
function satisfying 

r)C 

— -CG = 6{t-Hm-t>) (5) 



with G — > in |£ — £'| — > oo, where 

- a 2 



+ ^o(e))--?. (6) 



<9£ 2 vruvs " 4 
Formally, G reads 

G{tt\?J) = «o(0«o(0 +X)e- A " ( *- t ' ) «n(0<(0, (7) 

where £wo = 0, and £w n = \ n u n - The summation symbol, which may imply 
appropriate integration, is over the spectrum other than the point spectrum {0}. 
Because the system is translationally symmetric, u$ oc 

e CoC/2 V>o(0- Due to the 

known stability of the propagating wave front, the operator C is dissipative, so 
is the least upper bound of its spectrum. Hence, only uq contributes to the 
secular term (the term proportional to t — to) in Sip. Thus we can write 

61> = -(t -to) 6c1>' (Z) + (61>) r , (8) 

where (Sip) r is the bounded piece (regular part), and 

oc = — hm —. . (9 



This formula is easily justified if is isolated from the essential spectrum of L. 
In the pulled case (and for some examples of the pushed case), this condition is 
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not satisfied, ' but a more detailed argument shows that (9) remains valid. 
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One might immediately guess that 5c is the 0(e) change in the front speed, 
but the naive perturbation theory is not controlled, due to the secular divergence 
as to — > — oo. This divergence can be controlled by perturbatively renormalizing 
the constant of integration: xq = x R + Z, where x R is the finite observable 
counterpart of xq, the renormalization constant Z = a n(to, fi)e n and the 
coefficients a n are chosen order by order in e to eliminate the secular divergence. 
The quantity fi parameterizes the family of solutions to the unperturbed equation; 
it corresponds to the arbitrary renormalization point in the Gell-Mann — Low 
RG 24 

To order e the solution ip is given by 

i;(x, t) = MZo) - 8c(t - * W(£o) + (<Wr + 0(e 2 ) (10) 
= MO + £ anPo'iO ~ Sc(t - t )M(O + (^)r + 0(e 2 ), (11) 

where ^ = x — cot + x R (pt). Thus we can choose ea\ = (ji — to)5c to eliminate the 
divergence to 0(e). Requiring that ifj be independent of fx gives the RG equation 

-<**>. (-) 

Thus the speed of the renormalized wave is indeed cq + 5c. The formula (9) 
can also be obtained from the solvability condition for the first order correction 
5ip, and is an example of a very general relation between renormalizability and 

23 

solvability. Furthermore, (12) corresponds to the amplitude equation describing 

23 

the slow motion. This relation is also quite general. 

As an illustration of the use of the renormalized perturbation theory consider 
the following examples. The first, a pulled case example, is equation (1) with the 



nonlinear operator F = — ip) and the perturbation 5F = eip(l — ift). In this 

trivial case, the exact result is, of course, c* = 2\/T+~e, whereas (9) gives c* ~ 

2 + e. A more interesting pushed case example is provided by equation (3) with 

b E (0, 1/2). When 7 = 0, we have c*(0) = ^/^b+l/V^b. For non-zero 7, (9) gives 

c *( 7 ) = c *(0)-7c*(0) 3 s 4 (2s 2 + l)/10withs = 2c*(0)/(c*(0) + ^(O) 2 - 4). This 

agrees well with numerical calculations. For example, this result gives c*(0.08) ~ 

20 

2.696 for b = 0.1, while the corresponding value determined numerically is 
2.715. 

The perturbation theory result (9) can also be used to calculate heuristically 
the selected speed of the unperturbed system, using the structural stability idea. 
Within perturbation theory, a necessary and sufficient condition that c* be the 
selected speed is that 5c(c*) be bounded. For example, when F = — ip), the 
change in the velocity 8c(c) is zero as \\SF\ \ — > for all perturbations SF, which 
are both p-small and differentiable at the origin, only if c = c* = 2; for c > c* 
there exist such perturbations for which Sc does not vanish. A simple example 
of the latter is the perturbation SF = 9{u — A) (it — A)(l — it) — it(l — it), as 
A -> 0+. 

What is the physical significance of structural stability? Returning to the 
fuse analogy introduced above, we can imagine the fuse to be covered with a very 
thin film of water which quickly evaporates when heated to a temperature slightly 
above ip = 0; nevertheless the film suppresses tip ignition. Thus even a small 
perturbation can destroy (or drastically alter) the tip of a propagating front. For 
this reason, any front whose behavior is determined by its tip can be destroyed by 
such a perturbation. If and only if a front's behavior is independent of the details 
of its tip, can it survive such a perturbation and be structurally stable; hence, 
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a front with c > c* is not structurally stable, because no deviation is permitted 
from the required decay ahead of the front. On the other hand, the front with 
c = c* is insensitive to its tip, as we can see from our explicitly dynamical RG 
calculation. There, the leading edge is determined by the initial conditions, is 
not universal, and vanishes as t — > oo. Nevertheless, for sufficiently rapid leading 
edge decay in space, the asymptotic speed is c*: thus c* is independent of the 
details of the leading edge, so that the selected front is structurally stable. 
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